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SECTION  I 
INTRODUCTION 

The  purpose  of  the  program  STREAM  Is  to  theoretically  test  light  transmission 
through  complex  cyllndrlcally  symmetric  SHIVA  chamber  geometries.  (SHIVA  is  a 
Z-pInch  device  used  to  produce  hot  dense  radiating  plasmas  by  the  Ionization  and  J 
x B Implosion  of  a short  cylindrical  foil.)  If  sufficient  UV  light  Is  transmitted 
from  the  Ionized  foil  to  the  chamber  Insulator,  the  resulting  Ionization  of  the 
Insulator  may  cause  the  chamber  to  break  down  electrically. 

The  program  STREAM  has  the  following  characteristics: 

1.  Straightforward  input. 

2.  Three-dimensional  cylindrical  geometry  with  theta  symmetry. 

3.  Specular  or  diffuse  reflections  or  any  combination  of  the  two  phenomena 
may  be  specified.  Diffuse  reflections  can  be  treated  In  several  ways; 
for  example,  any  distribution  (Isotropic  or  cosine  distributions  are 
standard)  Is  allowed  and  It  can  be  oriented  around  the  surface  normal 
or  along  the  specular  reflection  angle. 

4.  Absorption  Is  calculated  for  each  reflection  and  It  is  a function  of 
the  angle  of  Incidence,  Index  of  refraction  and  extinction  coefficient. 
(Values  for  n and  k are  found  In  the  American  Institute  of  Physics  Hand- 
book for  various  wavelengths  and  materials.) 

5.  The  output  includes  relative  intensity  crossing  designated  surfaces  In 
the  geometry.  The  Input  Intensity  is  unity.  An  optional  output  file 
Is  available  for  2-D  microfiche  plotting  of  successful  (hits  on  the  in- 
sulator) photon  paths. 
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SECTION  II 
STREAM  INPUT 

The  Input  for  STREAM  Is  free  format  where  a comma  begins  a card  of  Input  and 
a slash  ends  It.  There  Is  a maximum  of  five  Input  sections.  The  first  two  are 
the  LAB  (for  label)  and  UNI  (for  units)  sections  which  are  each  very  short.  The 
LAB  section  Inputs  two  title  cards  which  may  be  changed  for  each  run  or  series  of 
runs.  The  UNI  section  Is  optional.  It  allows  Input  In  units  of  inches  (standard 
units  are  centimeters).  This  unit  conversion  applies  only  to  the  chamber  geom- 
etry card.  The  third  section  Is  the  SOU  (for  source)  section.  It  contains  the 

bulk  of  the  Input.  The  source  card  Is  followed  by  nine  or  ten  different  input 
cards  (a  card  In  this  context  contains  a stream  of  Input  terminated  with  a slash; 
several  IBM  cards  can  be  used  for  each  "card"  defined  here).  The  last  two  sec- 
tions are  the  END  and  STO  (for  stop)  sections.  Input  for  the  five  sections  is 

outlined  on  the  succeeding  pages.  Essentially  the  entire  problem  setup  Is  con- 
tained In  the  SOU  section.  (Note:  Only  the  first  three  characters  of  a header 

card  are  Interpreted;  the  rest  of  the  card  may  be  used  for  comments.) 

LABel  STREAM  Program 

,0,1/ 

first  title  card 
second  title  card 
UNIts 

,1,0, 0,1/  Chamber  geometry  is  Input  in  inches. 

SOUrce 

The  first  source  card  Inputs  several  integer  parameters  for  the  problem 
setup: 

,1,3,  NRMAX,  NWAVES,  NCCMAX,  NCSMAX,  NCRMAX,  NTIMES 

.NDMAX,  NPLOT,  MAXRAY,  NSTEPS,  NTOTAL , NORDER,  NSCAT,  NMMAX,  NTORUS/ 

The  above  variables  are  defined  as  follows: 
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NRMAX  — 


NWAVES  — 
NCCMAX  — 


NCSMAX  - 
NCRMAX  — 


NTIMES  - 
NDMAX  -- 
NPLOT  - 

MAXRAY  - 
NSTEPS  - 

NTOTAL  - 
NORDER  — 


NSCAT  - 


the  number  of  R-Z  cells  In  the  geometry  (each  cell  must  be 
trapezoidal  and  have  all  the  Interior  angles  less  than  180°). 

The  2-D  geometry  Is  rotated  around  the  Z-axis.  Possible  values 

^ 1. 

the  number  of  wavelengths  considered  In  each  property  set. 

Possible  values  1 ^ NWAVES  ^ 5 

the  number  of  property  sets.  Each  property  set's  results  are 
printed  across  a page.  Succeeding  property  sets  are  printed  on 
new  pages.  NWAVES  x NCCMAX  is  the  total  number  of  wavelengths 
and  materials  considered.  Additional  property  sets  do  not  in- 
crease the  running  time,  since  the  Monte  Carlo  is  done  once  and 
the  amount  of  absorption  per  reflection  is  the  only  difference 
between  different  wavelengths/materials.  Possible  values  ^1 

the  number  of  points  used  to  describe  the  source  polar  distri- 
bution. Possible  values  ^ 2 

the  number  of  points  used  to  describe  the  polar  distribution  of 
diffuse  reflections  relative  to  the  specular  reflection  angle  or 
the  surface  normal.  Possible  values  ^2 

the  number  of  time  bins 

number  of  detectors.  Possible  values  ^1 

logical  number  of  the  output  tape  used  for  creating  microfiche 
plots  (0  = no  plots).  NPLOT  = 9 implies  a TAPE9  file 

the  maximum  number  of  photon  paths  to  be  put  on  a plot  tape 

the  maximum  number  of  paths  through  which  each  particle  is 
followed  [(NSTEPS  - 1)  is  the  number  of  reflections  allowed.] 

the  total  number  of  histories  or  particles  followed 

the  number  of  "order  of  scatter"  edits.  If  NORDER  = (NSTEPS  - 1), 
then  all  possible  scoring  paths  are  edited  and  printed.  Possible 
values  ^ 0 

diffuse  distribution  reference  direction: 

0 = specular  reflection  direction,  1 = surface  normal 
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1 

NMMAX  — the  number  of  surface  materials  allowed  In  the  chamber  geometry. 

Possible  values  ^ 0.  If  NMMAX  = 0 or  1 , one  material  is  speci- 
fied. If  NMMAX  ^ 2,  then  additional  Input  Is  required  after  card 
5 (see  remark  before  card  6 explanation).  Also,  a card  after  card 
9 Is  used  to  describe  the  surfaces  containing  materials  2 through 
NMMAX  In  the  chamber  geometry  (see  card  10  remarks). 

NTORUS  --  this  card  Indicates  how  many  curved  surfaces  there  are  in  the 

geometry.  Each  curved  surface  Is  a portion  of  a circle  in  cross 
section  and  a portion  of  a torus  In  three  dimensions.  For  each 
curved  surface  Indicated  by  NTORUS,  special  Input  is  required  in 
card  9 (see  card  9 for  a description). 

The  second  source  card  describes  the  source  plane  and  the  source's  reference 
direction: 

,RMIN,  RMAX,  ZMIN , ZMAX,  RDIR,  O.ZDIR/ 

Definitions  of  these  variables  are  as  follows: 

RMIN  --  minimum  radius  of  source  plane. 

RMAX  --  maximum  radius  of  source  plane. 

ZMIN  — minimum  height  of  source  plane. 

ZMAX  --  maximum  height  of  source  plane 

RDIR  — R-axIs  direction  cosine  of  reference  direction. 

0 --  theta  direction  cosine 

ZDIR  --  Z-axIs  direction  cosine  of  reference  direction. 

NOTE:  RMIN,  ZMIN,  and  ZMAX  are  In  cm,  and  the  following  constraint  applies: 

[(RDIR)2  + (ZDIR)2]0,5  = 1.0 

The  third  source  card  describes  the  source's  relative  angular  distribution: 

ui  .Pi  .U2.P2 PNCSMAX  ,PNCSMAX/ 

is  the  cosine  of  the  angle  relative  to  the  cosine  direction  for  which  the 
probability  of  propagation  is  p^.  Two  special  cases  are  available.  ,0,0, 1,1/  is 
a cosine  distribution  (at  90°  the  probability  is  zero;  at  0°  the  probability  is  1). 

,0,1,1,!/  is  an  isotropic  distribution  (at  90°  the  probability  is  1.0;  at  0°  the 
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probability  Is  1.0).  Any  other  source  distribution  can  be  approximated  by  includ- 
ing several  (instead  of  two)  cosine  and  probability  pairs.  The  code  interpolates 
intermediate  values. 

The  fourth  and  fifth  source  cards  are  input  in  pairs.  They  are  repeated  for 
each  property  set  (e.g.,  if  there  are  three  property  sets,  three  pairs  of  cards  4 
and  5 must  be  included  in  the  deck.) 

Card  4 contains  the  mirror  reflection  fractions  for  NWAVES  wavelengths: 

,0.0,  0.0,  0.0,  0.0,  0.0/  mirror  fractions 

This  card  indicates  100%  diffuse  reflection  probability  for  each  of  five  wave- 
lengths in  a property  set. 

Card  5 contains  the  index  of  refraction,  n,  and  the  extinction  coefficient, 
k,  for  NWAVES  (five  in  the  following  case)  wavelengths: 

ni  ,kj  ,n2,k2,n3,k3,n4,k1|  ,n5,k5/ 

The  n and  k values  are  available  for  metals  in  the  American  Institute  of  Physics 
Handbook,  page  6-124.  If  NMMAX  is  2 or  greater,  the  above  sets  of  cards  4 and  5 
(mirror  fractions,  indices  of  refraction,  and  extinction  coefficients)  must  be  re- 
peated for  each  different  material. 

The  sixth  card  contains  the  diffuse  angular  distribution  which  can  be  calcu- 
lated relative  to  the  specular  reflection  direction  (NSCAT  = 0).  (Experimental 
evidence,  accumulated  with  sandblasted  copper  surfaces,  suggests  that  laser  light 
is  scattered  from  that  surface  ip  a cosine'squared-like  distribution  about  the 
specular  reflection  direction.) 

,Ul,Pl WNCRMAX’PNCRMAX/ 

The  meaning  of  y and  p is  the  same  as  on  source  card  3.  The  special  cosine  and 
isotropic  distributions  are  both  available. 

The  seventh  card  contains  the  detector  information.  The  detectors  are  de- 
fined on  any  boundary.  They  are  defined  by  cell  number  and  boundary  index  (1,  2, 
3,  or  4).  (The  cell  and  boundary  concepts  are  explained  later.) 

,NRMAX,4,1 ,1/ 

The  first  detector  is  on  the  fourth  boundary  cf  the  last  cell,  and  the  second  de- 
tector is  on  the  first  boundary  of  the  first  cell.  (NOTE:  Any  radiation  incident 

on  the  fourth  boundary  of  the  last  cell  is  considered  to  be  absorbed  in  the  insula 
tor. ) 
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The  eighth  card  of  the  source  section  contains  the  time  bin  information. 

There  should  be  one  more  time  boundary  than  NTIMES,  the  number  of  time  bins. 

,0.,  . . . time  bin  limits 

The  ninth  source  data  card  is  the  geometry  card.  The  geometry  consists  of  a 
series  of  adjacent  trapezoids  (none  of  which  Is  allowed  to  have  an  interior  angle 
greater  than  180°  or  overlap  adjacent  trapezoids).  It  is  input  by  pairs  of  (r,z) 
coordinates,  two  pairs  (called  quads)  at  a time.  Each  quad  describes  the  endpoints 
of  a line.  The  quads  are  input  in  order  from  one  end  of  the  geometry  to  the  other. 

Trapezoids  (called  cells)  are  automatically  formed  when  the  code  connects  the  top 
and  bottom  points  of  adjacent  quads. 

The  side  described  by  the  first  (in  order)  quad  is  side  1 of  the  cell.  The 
side  described  by  the  second  quad  is  side  4 of  the  cell.  The  bottom  border  of  the 
cell  is  side  2,  and  the  top  border  is  side  3.  There  must  be  one  more  quad  than 
cells  specified  by  NRMAX.  The  three-dimensional  characteristics  of  the  geometry 
occur  because  the  assemblage  of  trapezoids  is  symmetric  about  the  line  R = 0 
(theta  symmetry).  The  easiest  way  to  input  the  geometry  is  by  quads: 

,rlow’zlow’rMgh’zh1gh  (,st  face> 

,rlow’zlow’rhigh’zhigh  <2nd  facc> 


,rlow’z1ow’rhigh’zh1gh  (last  face'  NWWtl  quadl 

A diagram  of  the  trapezoidal  input  approach  is  in  figure  1.  The  cell  surfaces 
are  numbered  1 through  4. 

Figure  2 shows  a complete  geometry  with  one  baffle  inserted.  In  this  figure, 
particles  which  scatter  back  through  side  1,  cell  1,  are  lost.  Particles  which 
Impact  side  4 of  the  last  cell  are  successful  photon  strikes  on  the  insulator  and 
they  are  terminated.  The  associated  energy  is  tallied  for  each  wavelength. 


J 
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Figure  2.  Total  Geometry  (Example). 

An  additional  option  has  been  added  to  the  STREAM  program.  This  option 
allows  curved  surfaces  (arcs  of  circles  in  two  dimensions  or  portions  of  a toroid 
In  three  dimensions)  to  be  included  In  the  geometry  boundaries.  The  number  of 
curved  surfaces  to  be  Input  was  specified  by  NTORUS  on  the  first  SOUrce  card.  A 
curved  surface  is  Input  using  a special  quad  as  shown: 

.NFLAG,  NPOINT,  R',Z’ 

where  NFLAG  = 0 and  indicates  that  this  is  a special  quad. 


1 

1 
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NPOINT  = ±1,  and  +1  Indicates  an  arc  will  be  drawn  through  the  R',Z'  pair 
and  the  second  R,Z  pairs  of  the  preceding  and  succeeding  quads.  -1  indicates  the 
first  R,Z  pairs  of  the  preceding  and  succeeding  quads  will  be  used.  When  both 
sides  2 and  3 of  a cell  are  curved,  two  special  quads  need  to  be  inserted,  one 
with  NPOINT  positive  and  one  with  NPOINT  negative.  An  example  is  useful: 


,3.625,3.49,3.625,3.9  (Normal  quad) 

,0,+l  ,3.75,3.775  (Special  quad) 

,3.875,3.49,3.875,3.9  (Normal  quad) 

The  above  sequence  of  cards  indicates  that  an  arc  will  be  drawn  through  the  points 
(3.625,3.9),  (3.75,3.775),  and  (3.875,3.9). 

By  default  all  surfaces  in  the  geometry  have  the  optical  properties  of  mate- 
rial 1.  If  more  than  one  material  is  indicated  by  NMMAX  ^ 2,  then  specific  sur- 
faces may  be  selected  to  have  the  optical  properties  of  materials  2 through  NMMAX 
on  card  10. 

Card  10  follows  the  geometry  section.  It  is  an  optional  card  and  the  last 
card  in  the  SOUrce  section  of  input.  Input  is  in  pairs.  An  example  is  helpful: 

,-13,2  (±  cell  index,  material  index) 

,14,2/ 

The  cell  Index  indicates  the  cell  where  the  material  will  be  changed.  The  cell 
Index  is  preceded  by  a sign  which  flags  surface  2 or  surface  3.  If  it  is  negative, 
surface  2 Is  Indicated;  if  it  is  positive,  surface  3 is  indicated.  The  second 
value  Is  the  material  index  which  applies.  The  above  example  indicates  that  mate- 
rial 2 optical  properties  will  be  on  surface  2 of  cell  13  and  surface  2 of  cell  14. 

To  terminate  STREAM  input,  two  more  input  sections  are  used: 

END 

STOP 

,1/ 

The  following  guide  will  be  useful  to  users  of  STREAM. 
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STREAM  INPUT 

LAB  el 

,0,1/ 

first  title  card 
second  title  card 

UNIts 

,1 ,0,0,1/  for  geometry  section  (only)  In  Inches 

SOUrce 

,1,3 

,NRMAX,NWAVES,NCCMAX,NCSMAX,NCRMAX,NTIMES 
, NOMAX ,NPL0T .MAXRAY .NSTEPS .NTOTAL .NOROER .NSCAT .NMMAX , NTORUS/ 
,RMIN,RMAX,ZMIN,ZMAX,RDIR,0,ZDIR/  Source  plane  and  reference  direction 

,Ui,Pi,U2>P2>  • • •»yNCSMAX,pNCSMAX/s0urce  angu1ar  distribution 

,fi,f2,  • • •»fNWAVES/nv,rror  fract1on 
,n1,k1,n2,k2,  . . • •nNWAVES*kNWAVES/ 

Repeat  the  last  two  cards  If  NCCMAX  > 1 

Repeat  the  entire  combination  of  cards  4 and  5 If  NMMAX  ^ 2. 

,ui»P1»u2,P2,  • • • ’^NCRMAX’PNCRMAX^1^1156  reflect1on  angular  distribution 
.NRMAX, 4,1 ,1/two  detector  locations 

,0,  . . • ,tNTIMES+l/,time  bln 
C Geometry 

,I'low’zlow’rh1gh'zh1gh  (,st  face> 


’rlow’zlow’rh1gh’zh1gh  (,a!t  face> 

Insert  NTORUS  cards  In  the  above  geometry  with  the  form: 

*0,±,r' ,z 1 ; 

If  NMMAX  > 2,  Insert  card  with  as  many  pairs  as  required. 
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,±  cell  number,  material  number 


,±  cell  number,  material  number/ 
END 
STOP 
,1/ 


1 
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SECTION  III 

INTERPRETATION  OF  OUTPUT 

STREAM  produces  both  printed  and  plotted  output.  The  plotted  output  comes 
out  on  microfiche.  The  printed  output  Includes  the  geometry  outline  and  the 
printed  plots  of  several  sample  particle  tracks.  These  print  plots  are  In  pairs: 
r-z  coordinates,  then  r-0  coordinates.  A crossing  summary  which  describes  how 
much  energy  entered  and  exited  each  surface  of  each  cell  In  the  geometry  is  in- 
cluded. The  last  set  of  output  Is  the  most  Important;  It  has  detector  results 
listed  both  by  time  bln  and  by  number  of  collisions  needed  to  score  at  that  de- 
tector. This  output  Is  In  column  pairs.  The  left  column  of  each  pair  has  the 
outputted  result  for  each  wavelength,  and  the  right-hand  column  has  the  corres- 
ponding error  bar.  At  the  top  of  each  pair  of  columns,  the  mirror  fraction  Is 
listed  along  with  the  Index  of  refraction  and  extinction  coefficient  for  the  wave- 
length and  material.  Totals  are  tabulated  for  each  column.  The  order  of  the  "By 
Collision"  output  Is  determined  chronologically  in  the  code  by  the  order  in  which 
successful  hits  on  each  detector  occur.  Figure  3 has  a sample  of  the  order  of 
scatter  edit.  Figure  4 Is  a sample  microfiche  plot  of  a z-cross  section  chamber 
with  three  baffles. 
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STREAM  PROGRAM 
MASTER  FILE  < 75) 


FOUR  BAFFLES— STRAIGHT  CHAMBER 
100*  SPECULAR  DISTRIBUTION 


OUTPUT  FOR  CELL  AND  SIDE  21 


4. TIME  TRANSLATES  2.6121E-10 


TIME  dependent  SET 


IIN 

Time 

1 =W  AVE 

ERROR 

?=WAVE 

ERROR 

l.OCE*Ofl= 

MIRROR 

1 »00E»0O= 

MIRROR 

l.P5E»00= 

N 

9, 50E-0 1 = 

N 

7.00E-0!= 

K 

7.80E-01= 

K 

1 

1.00E-99 

7.73E-C9 

1.63E-C1 

4.18E-08 

1.57E-G1 

2 

2.00E-99 

6.75E-22 

4.47E-01 

1.39E-19 

4.11E-01 

3 

3.00E-09 

0. 

0. 

0. 

0. 

TOTALS 

7.73E-09 

1.63E-01 

4.18E-08 

1.57E-01 

1 0Y 

COLLISION 

13 

U68E-13 

1.49E-01 

3.79E-12 

1.48E-01 

19 

9.25E-18 

1.37E-01 

3.74E-16 

1.32E-01 

15 

3.73E-15 

1.69E-01 

1.33F-13 

1.65E-01 

10 

3.61E-17 

1.30E-01 

2.35E-15 

1.27E-01 

16 

2.60E-15 

9.22E-02 

1.08E-13 

9.03E-02 

12 

1.82E-12 

1.41E-01 

3.05E-11 

1.39E-P1 

29 

7.08E-27 

3.87E-01 

5.82E-24 

3.32E-01 

21 

4.61E-20 

1.99E-01 

6.7IE-18 

1.91E-01 

24 

6.18E-22 

2.09F-01 

1 .47E-19 

1.94E-01 

9 

4.65E-10 

1.31E-01 

3.84E-09 

1.31E-01 

14 

l.l^E-13 

9.08E-02 

2.95E-12 

8.95E-02 

22 

1.10E-2J 

2.94E-3 1 

1 .68E-18 

2.72E-01 

23 

1.50E-21 

2.56E-71 

3.21E-19 

2.38E-01 

28 

1 « 1 IE-25 

6.31E-01 

6.C6E-23 

5.75E-01 

7 

7.12E-09 

1.77E—31 

3.70E-08 

1.77E-01 

11 

1 .48£- 1 2 

2.37E-  U 

2.12E-11 

2.36E-01 

17 

1.88E-16 

1.38E-M 

KOlE-14 

1.36E-01 

25 

2.42E-23 

3.32E-0 1 

8. 16E-21 

3.1 3E-9 1 

30 

1.I2E-26 

5 »66r—C l 

7.R7E-24 

5.24E-01 

26 

4.83E-24 

4.25E-C  1 

1.77E-21 

3.90E-C1 

20 

5.20E-19 

2.16E-M 

5.15E-17 

2.08E-01 

27 

3.75E-25 

3.87E-01 

2.1  IE-22 

3.59E-01 

8 

1.4CE-1G 

5.87E-01 

9.20E-10 

5.73K-01 

10 

6.69E-13 

5.77E-01 

8.07E-12 

5.77E-01 

TOTALS 

7 . 73E-09 

1.63E-01 

4. 18E-08 

1.57E-01 

Figure  3.  A Sample  of  a Typical  Order  of  Scatter  Edit. 
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SECTION  IV 
STREAM  ACCESS 


The  program  STREAM  may  be  accessed  three  ways.  A permanent  file  of  the  STREAM 
OLDPL  Is  available  on  the  CDC  7600  computer,  a deck  of  cards  Is  available  in  the 
X-Ray  Simulation  Branch's  terminal  room,  and  a copy  of  the  STREAM  OLDPL  Is  avail- 
able on  tape  BF98  which  has  the  label  STREAM. 

To  access  STREAM  and  run  It,  the  following  control  cards  are  used: 

JOB  CARD 
ACCOUNT  CARD 

ATTACH  (OLDPL, STREAM, ID=DYSXTMJ) 

UPDATE  (K,P=0LDPL,L=0) 

REWIND  (COMPILE) 

ATTACH  (FTNC0MP4P2L383) 

LIBRARY  (FTNCOMP) 

RFL  (160000) 

FTN  (A,I=C0MPILE,0PT=2 ,L=0) 

REDUCE. 

LIBRARY. 

REQUEST  (TAPE9  ,*PF) 

LDSET  (MAP=0) 

LGO. 

CATALOG  (TAPE9 .CHAMBER , ID=DYSXJAJ  ,RP=3) 

7/8/9 

* I DENT  FIXIT 

♦COMPILE  MAIN, BLOCK, BLOCKV,ARGSET,GETARG,DUPEMA,DUPEZA,FILLMA,FILLZA 
♦COMPILE  SWAPMA .SWAPZA ,ZIN ,ZLOG ,ZOS ,ZQRT,ZTAN2,ZXP .MOVEMA .ZEROUT 
♦COMPILE  CARDIN .DIGIT4 .DUMPIT .ERRORS , INTVAL , INTWRD , IOLINE .KOMMON .LABEL 
♦COMPILE  LABOUT,LASTEQ,LIMVAL,LINSET,MOVETT,PACKUN,RANNO,ROTATC, ROTATE 
♦COMPILE  ROTATX, SCALER, VECTOR, DEFINE, MASTER, READIT, HEADER, IARRAY .ICHECK 
♦COMPILE  JCHECK,KCHECK,TABCON,TOCSET, INDEX, LABSET,BOOTIT,DATBAS,DATRUN 
♦COMPILE  BYLINE, STREAM, HELIX 
7/8/9 

INPUT  DECK 
6/7/ 8/9 

The  tape  9 file,  which  Is  cataloged  at  the  end  of  the  control  card  STREAM, 
contains  the  raw  data  which  are  used  for  generating  two  dimensional  microfiche  plots. 
The  tape  9 file  was  specified  by  NPLOT  on  the  SOUrce  option  card. 
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SECTION  V 
MICROFICHE  PLOTS 


A program  called  FISHY  Is  used  to  generate  the  microfiche  plots  of  photon 
paths  In  the  chamber  geometry.  This  program  gets  Its  Input  from  the  permanent  file 
CHAMBER  which  contains  the  TAPE9  file  created  by  STREAM.  It  needs  the  DYSPLOTLIB 
to  use  in  generating  the  microfiche  plots. 

The  following  control  cards  are  needed  to  run  FISHY: 

JOB  CARD 
ACCOUNT  CARD 
FTN(A,0PT=2,L=0) 

ATTACH (DYSPLOTLIB  ,ID=DYSXBBB ) 

LIBRARY  (DYSPLOT) 

ATTACH  (TAPE! , CHAMBER,  ID=DYSXJAJ  ,CY=1 ) 

REWIND  (TAPE1) 

FILE  (JAOFIS,RT=S) 

LABEL  (JAJFIS,L=CHAMBER) 

REQUEST  (JAJFIS ,PE,E,VSN=HG71 ) 

MAP  (ON) 

LGO. 

7/8/9 

FISHY  DECK 
6/7/8/9 

In  the  above  case,  the  microfiche  tape  that  Is  being  generated  Is  HG71 . Any 
nine-track  tape  may  be  substituted  on  the  REQUEST  card. 

A listing  of  FISHY  follows: 
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